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Abstract 

We present a definition of angular momentum for radiative spacetimes which 
does not suffer from any ambiguity of supertranslations. We succeed in providing 
an appropriate notion of intrinsic angular momentum; and at the same time a 
definition of center of mass frame at future null infinity. We use the center of mass 
frame to present the asymptotic structure equations for vacuum spacetimes. 

PACS numbers: 04.20.-q, 04.20.Cv, 04.20.Ha 

1 Introduction 

The difficulty associated with a sensible definition of angular momentum at future null 
infinity T + (scri) is related to the complicated structure of the Bondi-Metzner-Sachs (BMS) 
group of asymptotic symmetries. There is a multiplicity of Poincare subgroups, one for 
each supertranslation. 

Because of this difficulty one can find in the literature several inequivalent definitions of 
angular momentum |TW66] |Win 68j [Br a75] |Pri77j |Str78j [ Pen82j |DS84j which suffer from 
the so called supertranslation ambiguities. The only definition of angular momentum 
without supertranslation ambiguities is the one we presented in reference [Mor86j; however 
our construction can be criticized on the grounds that we select a unique frame based on 
the properties of future null infinity in the limit for the retarded time going to — oo. A 
global choice of this nature can hardly be thought to be able to describe the local meaning 
of intrinsic angular momentum of the sources at an arbitrary retarded time of a radiating 
system. 

In order to provide with an unambiguous notion of angular momentum in radiative 
spacetimes, it is essential to make use of a reference frame system, which embodies the 
notion of rest frame [WinSOj . 
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The problem of the characterization of the notion of a rest frame for isolated systems 
has been the matter of concerned for many years. Early attempts used the properties 
of an optical parameter of a congruence of null geodesies reaching future null infinity; 
namely, the shear. When the spacetime is stationary, it is known |RP68j that one can 
choose a Bondi system for which the shear <r is zero; where the freedom in this choice is 
just a translation. The sections for which the shear is zero are known as 'good cuts'. We 
see then that for stationary spacetimes there is a simple geometrical way to get rid off 
the supertranslation problem. However for radiating spacetimes things are not so easy. 
In presence of gravitational radiation, the good cut equation has in general no solution; 
but it was suggested in the literature that the solution to this can be obtained if X + 
is complexified. This approach has lead to the beautiful construction of TC spaces(See 
reference |NT80 and references there in.). Due to the fact that Ti spaces are intrinsically 
complex, it remained the difficulty of the interpretation, in this construction. 

This situation motivated us in the past to introduce the so called 'nice sections' [ Mor88| 
at X + . These sections are defined in terms of the nice section equation; which recently 
has been proved MD98J [DMOOJ to have solutions in terms of a 4-parameter family of 
translations. Furthermore the solutions have the expected physical properties of a rest 
frame[DM00j; for example, having a nice section So, all other nice sections Sf obtained 
from future timelike translations happen to be to the future of 5*0. 

The nice section construction singles out precisely, in an intrinsic way, a Poincare 
structure from the infinite dimensional BMS group. In particular, given a fixed obser- 
vational point p at X + , there is precisely a 3-degree of freedom of spacelike translations 
which generate all the nice sections that contain p. In contrast, without this construc- 
tions there is an infinite dimensional family of general sections that contain p, one for 
each supertranslation. 

The success of the nice section construction allows us to define a well behaved notion 
of intrinsic angular momentum along with the associated notion of center of mass frame. 

In section El we briefly review the definition of nice sections. We make use of charge 
integrals of the Riemann tensor to present in section |3] both the angular momentum 
definition and the center of mass frame. In section H] we study the asymptotic structure 
equations of a general radiating spacetime in terms of the center of mass frame. We briefly 
compare our presentation with previous definitions of angular momentum and make final 
comments in section 03 

We will make use of the OHP |GHP73] notation. 

2 Rest frame systems at future null infinity 

If a physical system can be ascribed a Poincare structure, one should be able to define its 
momentum vector P a and angular momentum tensor J ab ; and then one would define a rest 
frame by demanding that the total momentum has no spacial components in this frame. 
However an isolated gravitating system has the BMS structure [Mor8 7j; therefore it is 
natural to define a rest frame by transforming to a section for which the supermomentum 
has no spacial components. This brings us to the issue of which supermomentum to use. 
In the nice section construction we use Mor88j the supermomentum, that for simplicity we 
are going to call 'psi'. The supermomentum psi has some interesting properties, but it is 
important to emphasize that it does not coincide with the Geroch supermomentum |( ler77j 



2 



or with the Winicour one[TWfifij [ Wrn68] jWTn80] |G W81 j . 

Given an arbitrary section S of Z + , one can, without loss of generality, choose a 
Bondi coordinate system (u, (,(), such that u = determines the section S. Then, the 
supermomentum psi on S is given by the expression 

p lm (S) = --L / yj m (c, CMu = o, c, C)^ 2 , (i) 



where (iS* 2 is the surface element of the unit sphere on S, Yi m are the spherical harmonics, 
the scalar \1/ is given by 

■q? = + + g2- 0; (2) 

where ^ * s ^ ne leading order asymptotic behavior of the second Weyl tensor component, 
do is the leading order of the Bondi shear respectively, where we are using S to denote the 
GHP[GHP73 edth operator of the unit sphere, and where a dot means partial derivative 
with respect to the retarded time u. 

The first four components of the supermomentum, namely the case I = and the 
three cases 1 = 1, determine the Bondi energy-momentum vector. More specifically 

(n = (Poo, -^( p n - *Vi)> ^( p n + Pi,-i), ^io) , (3) 
where a = 0, 1, 2, 3. 

All other nice sections can be obtained from S by the appropriate supertranslations 7. 



Therefore one can impose the nice section condition on 7. It was shown [Mor88 _MD 
that this condition can be cast in the following equation 

5 2 B 2 7 = *( 7 , C, C) + ^ 3 (7, C, C)M( 7 ), (4) 
where the conformal factor K can be related to the Bondi momentum by 

with 

(n = (1 1+1 c-c cc-i \ (6) 

1 ) V 'i + CC'*(i + CC)'i + CC/ 

and P a is evaluated at the section u = 7; which is calculated through the integral @. 
The rest mass M at the section is given by 

M = (7) 

where the indices are raised and lowered with the Lorentzian flat metric rj a b at scri as 
described in reference |Mor86j . 

As it was stated in the introduction, there exists a four parameter family of solutions 
of equation (jlj) with the expected physical properties MD98J [DM00 . One can always 
express 7(C ; C) i n terms of the translation part 7/(C,C) an d its proper supertranslation 
part 7n-(C,C)j where using spherical harmonics 7/ is expressed in terms of loo an d Y\ m \ 
while 7// is expressed in terms of Yi 2Tn with / 2 > 2. The result is that for each translation 
77 there is a proper supertranslation 77/. 
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We have assumed that the original section S determined by u = is a nice section. 
In other words 7 = is already a solution of (HJ), and so K(0,(,() = 1- For any other 
nice section 7 7^ in general the Bondi momentum will have spacelike components with 
respect to the original system, and so one would have K ^ 1. In that case, in order 
to observe that the space like components of the supermomentum are zero, one needs 
to make a Lorentz transformation of the BMS group, that aligns the generator of time 
translations with the Bondi momentum, at the new section. 

The need of supertranslations and Lorentz transformations for grasping the idea of 
rest frames, is what makes the discussion of radiative spacetimes so complicated; and it 
is clear that one single Bondi system does not provide with a rest frame unless the metric 
is stationary. 



3 Intrinsic angular momentum and center of mass 
3.1 Charge integrals of the Riemann tensor 

The rest frames provided by the nice sections construction permit us to define in an 
unambiguous way the physical quantities of an isolated system, in particular its angular 
momentum. 

It is convenient to approach the notion of physical quantities by the method of charge 
integrals of the Riemann tensor, as it was used in our previous work |Mo r86j. 

Given a 2-sphere S a charge integral of the Riemann tensor is a quantity of the form: 

Qs = [c (8) 

where the 2-form C ab is given in terms of the Riemann tensor by 

C ab = R* ab cd w cd} (9) 

and where the 2-form w a b is to be chosen appropriately. 

We will assume that the sphere S is actually a section of X + . Since the spacetime is 
asymptotically flat, it is possible to express the metric, in a neighborhood of scri, around 
a flat background, namely 

g a b = Vab + h ab , (10) 
where at scri, the conformal regular metric g ab is given by 

g ab \ I+ = £l 2 rj ab ; (11) 

in other words, at scri one has Q 2 h ab = . 

Let us consider, for a moment, the situation in which the expansion (fTT)|) can be 
extended to the interior of the spacetime; which, to simplify the discussion, will be as- 
sumed to contain no singularities. Let E be a spacelike hypersurface in the interior of the 
spacetime but that reaches asymptotically future null infinity; in such a way that in the 
conformaly completed spacetime, £ can be extended to scri with boundary 5'. Then, it 
is deduced that the charge integral on S can be expressed as an integral on E, by the use 
of Stokes' theorem, namely 

Qs= [ C= [ dC. (12) 
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(14) 



The exterior derivative of C can be expressed by 

dC abc = ^e abcd *R* def9 V e w fg ; (13) 

where it is interesting to note that a property of the double dual of the Riemann tensor is 
that its trace gives the Einstein tensor, namely: *R* abcd g bd — G ac = R ac — \g ac R. Therefore 
the previous equation can be written 

dC abc =^e abcd *R* def9 ^T e fg + ^9efV g - T^gegVfj 

= \e abcd {-2G d9 v g + *R* def9 T efg ) , 

where T abc is the traceless part of V a w bc and v c its trace; in other words 

S/ a w bc = T abc + ^g ab v c - ^9acV b ; (15) 

where it can be checked that From equation ((H)) one observes that if 

the vector v a where a Killing vector of the metric i] ab and T abc where 0(h), then the 
charge integral will give precisely the conserved quantities in the context of linearized 
gravity. This analysis ensures that this charge integrals admit the appropriate physical 
interpretations in the linearized gravity regime. 

In what follows, we will only assume that the spacetime is asymptotically flat at future 
null infinity [Mor 8 7j . 

Another interesting property of the double dual of the Riemann tensor is the one 
associated with the Bianchi identities, namely *R* d ^ e ^3\ = o, from which one can prove 
the relations [Mor86j 

T efg = ~ *R*^» (T {ef)g - T {eg)f ) ; (16) 

where it is important to note that the factor involving the tensor T abc have a very simple 
form when expressed in terms of the spinorial notation, that is 

2 

3 ( T (ef)9 - T (eg)f) = ^E'{EWfG) ^F'G' + c.c. ; (17) 

where c.c. means complex conjugate. 

Observing equation (pPfj) and recalling the analysis in the context of linearized gravity, 
it is natural to study at scri the conditions on w given by the equations 

-V/' w AB + c.c. = v BB ' (18) 

and 

Ve'(e wfg) = 0; (19) 

where the vector v BB ' is a generator of asymptotic symmetries. 

The generators of Lorentz rotations are a special case of asymptotic symmetries. The 
vector field v a , defined in a neighborhood of scri is said to be an asymptotic symmetry if 
it satisfies 

V ( .t%) = ^f, (20) 
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where S a b has a regular extension to I + . 

In general an asymptotic symmetry v a can be expressed by its components, in terms 
of a null tetrad frame 

v a = Vn t - Vmm a -v m rh a + v e n a ; (21) 

where we are using the standard complex null tetrad that satisfies i a n a = —m a fh a = 1 
with all other scalar products giving zero. As usual the tetrad vector £ a is chosen to 
generate a congruence of null geodesies reaching scri. A regular tetrad at future null 
infinity [l a , m a , fh a , h a ) can be constructed from the following relations 

t = Q 2 t, (22) 

m a = ttm a , (23) 

n a = h a ; (24) 

where Q can be taken as the inverse of the affine parameter r of the null geodesies reaching 
scri. 

Since the asymptotic symmetries are tangent to Z + , the tetrad components have the 
following behavior 

Vn = rv° n + v 1 n + 0(-), (25) 
r 

v m = rv° m + v 1 m + 0(^), (26) 



t,, = «0 + i + O(I). (27) 



The leading order behavior of the asymptotic symmetries is given by 

v° m = da, (2? 



v 



o 



1 /- B N 1 



„ 2 (pVn + 5v m ) = -ti8(a + a), (29) 

v? = x(C,0-u^m(a + a) (30) 

where we are using a Bondi tetrad Mor86 , x an d a are functions on the sphere with spin 
weight and satisfying: x = x, X = 0> a = and S 2 a = 0. Let us note that the function 
"a" only appears under the action of the operator edth; then, without loss of generality 
we can assume that it satisfies 55a = —a; which says that a can be expressed in terms of 
spherical harmonicas Y\ m with 1 = 1. 

The relation (j!8)) at scri can be expressed in terms of the spinorial components of the 
regular dyad 

w AB = w t A t B - wi (6 A l B + l a 6 b ) + w 2 o A o B , (31) 

by 

1 . . 

U>2 = ~^ v ^ ( 32 ) 
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1 . . 

w 1 + w 1 = -g^, (33) 

wx + wx = ~ (dw 2 + Qw 2 ) ; (34) 
while condition (fTTj|) at scri becomes 

dw 2 = 0, (35) 

w 2 = 0, (36) 

dw = —2 <7o w ii (37) 

ibx = ~dw 2 , (38) 

-w + Bwx + (To w 2 = 0. (39) 
It is interesting to see whether the solutions of (fTHJ) at Z + , namely (p?2~|l - (p5lj) also satisfy 

Before studying this for the stationary and radiating case separately, let us note that 
for any 2-form w defined on a sphere S, one can define the charge integral (JEJ) obtaining 

Qs(w) = 4 J [-w 2 {9fi - $io) + 2wx(^ 2 - $ n - A) - w {^ 3 - $21)] c^ 2 + c.c. (40) 

where (wq, Wx, W2) are the components of wab with respect to the dyad adapted to S, 
namely for which the vectors m and rh are tangent to S and £ is outgoing. If w has a 
regular extension to scri we can take the limit of S to a section of X + obtaining 



QsW 



4 J [-w 2 ^x + 2wx^° 2 - w ^l] dS 2 + c.c. (41) 



where the upper-script denotes the leading order behavior of the respective Weyl com- 
ponents. The charge integral at X + turns out to be independent of the Ricci tensor since 
as we have shown in [Mor87 j for any asymptotically flat spacetime, the components of 
the Ricci tensor go to zero faster than the accompanying Weyl components in the above 
terms. Furthermore, using that ^3 = — <3ao an d equation (}3*Tj) one can rewrite this as 



Q s (w) = 4 J [-w^ + 2w 1 (^ 2 + a a )] dS 2 + c.c. 



(42) 



It is important to note that the solutions of equations ([32 )1 - ([34 )1 can always be made to 
satisfy equations (|33|) - (p?%|) and so one can obtain expression for the charge integral, 
which is independent of the validity of equation (jHHJ)- 
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3.2 Stationary spacetime case 

For the case of stationary spacetimes one can solve the set of equations (J35j) - (J39|) with 
solution 

W2 = -\®a, (43) 

w 1 = w 1 °((,Cvo) + ]:uma, (44) 
o 

wo = w° °(C, C, °o) + u (-23< + ^ 8o) - X -u 2 3 2 8a (45) 

where 9 2 a = 0, u^ and are spin weight and 1 functions respectively that solve the 
equations 

1 - 1 3 
9 2 ^i° - o^°"o Sa + 7^0 3Sa = -5<t u> 2 - -cr <3u> 2 (46) 

and 

S^ = -2(70^?°. (47) 
Let us note that if one uses the potential 8 of the shear satisfying 

<J = <5% (48) 

then, the component W\ can be expressed by 

Wl = b+-dSda + -(u-5ma; (49) 
3 6 

where the spin quantity b satisfies b = and <5 2 b = 0. 

This procedure provides with a two-form with the functional dependence 



w 



o 

AB 



w °ab C,^o(C,C), a, b) . (50) 



Let us observe that from one of the Bianchi identities at scri, it is obtained 

m° 2 = ^°-2^°a , (51) 

but for a stationary spacetime ^5 = and also f || = in a Bondi system; therefore, 
one deduces that ^ is a constant. Then, since the nice section condition requires \1/ = 
^2 + oo^o + S 2 o"o = ^2 + 9 2 ^o to be a constant; one concludes that S 2 a = and therefore 
that (To = 0. In other words the nice section condition coincides in the stationary case 
with the good cut equation condition. However, as we have noticed, the nice section 
equation has solutions in the radiating case, while the good cut equation does not. 

In the stationary case we have that the charge integrals of the Riemann tensor at scri 
are given by 

Q S (w) = J C = A J {-w 2 + 2 Wl °) dS 2 + c.c. ; (52) 
where w 2 is given by (|33J) and W\ is 

Wi — b H — ttSSa. (53) 
6 
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Let us note that a involves 6 real constants associated with the Lorentz rotations, and 
that since ^ is a real quantity, b contributes to the charge integral Qs(w) with other four 
real constants associated with translations. 

We see then that this construction provides with the appropriate rest frames for sta- 
tionary spacetimes; and that the vector v a (a, b + b) calculated from the divergence of the 
2-form w, equation (|18j) . is a generator of BMS symmetries. 

The first term in the integrand of equation ()52|) includes the Weyl component 
which is known to describe the angular momentum in the Kerr geometry for example. 
In the second term we recognize the component ^ which form the supermomentum psi 
for this particular stationary case. Therefore equation (|52"|) has the same structure as 
the expression that, in special relativity, relates the angular momentum J ab , the intrinsic 
angular momentum S ah and the linear momentum P a , namely 

jab = £ab + fl apb_ pa R b ^ ^ 

where a, b are numeric spacetime indices. Given a rest reference frame in Minkowski 
spacetime one needs to use the spacelike translation freedom (R h ) appearing in expres- 
sion (p>f|) in order to single out the center of mass reference frame. In the center of mass 
frame one has J ab = 5* ab ; that is the total momentum coincide with the intrinsic angular 
momentum. Since the intrinsic angular momentum satisfies S* b Pb = 0, one can charac- 
terize the center of mass frame as that rest frame for which J 0i = 0; where i = 1, 2, 3 is a 
spacelike numeric index. 

In the nice section construction we have seen that given an observation point p at 
scri, there is a 3-family of rest frames, associated with the choice of spacelike translations. 
We also observe that the charge integral (|52"|) has the appropriate angular momentum 
behavior. Therefore we also need a condition analogous to the Minkowskian case to single 
out the center of mass reference frames. 

It can be seen that the condition to be imposed in the charge integral case is 

Qs( a ) — for all a = a; (55) 

where it is understood that one takes b = in this equation. The quantity a is in principle 
complex; so this condition makes use precisely of a 3-degree of freedom, which is associated 
with spacelike translations. 

This is the appropriate condition that leaves a one-dimensional family of nice sections 
S that can legitimately be called center of mass frames. We can see that the center of 
mass frames are generated by time translations in the nice section construction. 

Using these frames S cn , the intrinsic angular momentum j is defined through 

j(w) = Qs c »; (56) 

where to determine w one chooses a = —a and b = 0. 

The same charge integral can also be used to calculate the Bondi momentum, which 
requires to take a = and 6^0. 

This prescription singles out the center of mass frame for stationary spacetimes and a 
Poincare subgroup of BMS generators. 
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3.3 Radiating spacetime case 

In the presence of gravitational radiation, the shear ctq depends on the time u and the 
set of equations (|32|) - (}39|) in general have no solution. But it is interesting to note that 
one can define a 2-form w at scri, which in some sense is the closest to the stationary 
prescription, as we describe next. 

Let (u, C-C) be some Bondi coordinate system at future null infinity. A generator 
of scri can be identified with the line determined by £ = constant (£ = constant). 
Let us consider a point along a particular generator of scri, denoted by p(r), with r a 
monotonically increasing time parameter. 

Let us recall [DMOOJ that the set of nice sections form a four parameter (T, R) family 
that we now label S^ TI ^y, where {T,R) can be identified with a translation of the BMS 
group. 

Then, for a given fixed r, one has a 3-parameter (R) family of nice sections S^ T ^ 
that contain the point p(r); where all this family share the same T(r). Given one of this 
nice sections S, T g,, we can always identify it with the condition u = jg, where jsiCX) 
is the supertranslation that defines the corresponding section. On S, T fo we define the 
2-form wg n (u = 75, C> C> a i b) as the solution on S^ T ^ of equations (|T%|) and (fTTIjl for the 
radiation data as - (u,(,Q = a (js, (,()■ Since by definition as - does not depend 
on the time coordinate u, the 2-form wq ^ is the restriction on S trr of the solution 

' (T,R) (T,R) 

of the stationary problem when the radiation data is a$ T ~- More concretely, one can 
express the 2-form ws {Tji) in terms of its spinor components by the prescription stated 
in equations PHjl - (jTfjl . if we chose the Bondi system such that S^ T ^ is characterized by 
u = constant. 

In order to single out the center of mass section S cm from the 3-parameter (R) family 
of nice sections that contain the point p{r), we demand 

^ S {T,R) ^ = fOT fl = ( 57 ) 

This is the analog in the Minkowskian case of the condition J° l = 0; that we mentioned 
earlier. Let us note that the quantity a involves six real parameters; so that the condition 
a = a involves three degrees of freedom which are taking care of by the corresponding 
freedom involved in (R). In this way, for each r we select the unique S cm (T) which contains 
the observation point p(r) at scri. 

Given the radiation data o"o(w,C ; C)) this construction provides with a smooth DMOOJ 
one-parameter family of non-crossing center of mass sections S cm {r) along with a smooth 
2-form 

w(t, C, C; a, b) = w St (u = 7 5t , C, C; a > b) (58) 

on scri. 

The 2-form wab as defined by (}58|) will not satisfy all equations (J32|) - (}39|) : more con- 
cretely, most of the equations involving time derivatives will not be satisfied in general, 
but the error will be 0(&o). That is, if do = on a particular center of mass section S cm , 
then all equations (fTHj) and (fT^j) . or equivalently (p^ - ljBTIj) . are satisfied on S cm . 

Using these center of mass frames S cm , the intrinsic angular momentum j is defined 
through 

j(w) = QsM; (59) 
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where as before, in order to pick up the intrinsic angular momentum, one takes a = —a 
and 6 = 0. 

With this definition the charge integral of the Riemann on S cm , can be expressed by 



Qs cm (w) = f C = 4 I (-W2{*1 + 2a cV + 9 (a o a )) 

+ 2 Wl (*§ + o- <To + a 2 ab)) dS 2 + c.c. 

where we have used that 

(wi 5V ) dS 2 = I ( a Scfo + ^c3 (cr (f ) J w 2 dS 2 . 



(60) 



(61) 



It should be emphasized that in this way, given an asymptotically flat spacetime, 
one constructs a unique regular family of non intersecting sections |DM00j at future null 
infinity, which are the center of mass frames, and that can be used to describe the detailed 
multipolar asymptotic structure of the spacetime. 

In order to see whether our definition of angular momentum has the appropriate 
behavior in the presence of gravitational radiation, let us imagine a system in which one 
can distinguish three stages; starting with a stationary regime, passing through a radiating 
stage and ending in a stationary regime as depicted in figure [T] By construction it is clear 
that our definition gives the correct intrinsic angular momentum in the first and third 
stage; even though the center of mass frames of the first and third stage are related in 
general by a supertranslation. 

This is the only definition of angular momentum that satisfies this expected property. 



stationary stage 



radiating stage 




stationary stage 



Figure 1: Behavior of the notion of intrinsic angular momentum in a spacetime with three 
stages; beginning with a stationary stage, continuing with a radiating stage and ending 
with a stationary stage. 
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4 Asymptotic structure of asymptotically flat space- 
times 



4.1 Basic variables 

Having obtained the center of mass frames, we now take the opportunity to study the 
asymptotic structure in terms of a null tetrad which is adapted to the center of mass 
frames. 

We start by defining the basic variables that we will use. 

One can express the asymptotic geometry in terms of a complex null tetrad (£ a , m a ,m a , n a ) 
with the properties: 

g ab t n b = -g ab m a m b = 1 (62) 
and all other possible scalar products being zero, the metric can be expressed by 

9ab = L n b + n a £ b -m a m b - m a m b . (63) 

Newman and Penrose [NP62j showed that such a null tetrad is easily related to a dyad 
of spinors (o a ,l a ), and that they are very useful for the discussion of the asymptotic 
structure. 

Using the null polar coordinate system (x°, x 1 , x 2 , x 3 ) = (u, r, (C + £), |(£ — £)) one 
can express the null tetrad as: 

L = (du) a (64) 



dr 



(65) 



= (66) 

* =*(■&)' (67) 



with i — 2,3 and where the components £\ U and X % are: 



^=f+o(^), f 3 = | + 0(^), (69) 



with 

H = V2P V, Co = (70) 

(i+CC) 

sphere; 

Ui , „ / 1 



where V = V(u, (, () and the square of P = - ^ 1S the conformal factor of the unit 



U = rU 00 + U + — + O[-) 1 (71) 

where 



TJ V TJ l K TJ *° + *° (70) 
(7oo = y, Uq — — - Ay, (7i = , (72) 
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where Ky is the curvature of the 2-metric 

dS 2 = — Lj d( d(- (73) 

where the regular conformal metric restricted to scri is precisely g |j+= —dS 2 . In terms 
of the edth operator cV of the sphere (|73*j) the curvature Ky is given by 

K y = 1 g y 5 y y-^^ B ^ + ( 74 ) 

Finally, the other components of the vector n a have the asymptotic form 

* 2 = O (p) , * 3 = O (1) . (75) 

4.2 Asymptotic behavior of the spin coefficients 

The optical spin coefficients of the outgoing congruence of null geodesies are 

p = -;-^ + °(1) (*» 

and 



ty £i ly'-i ' qy 

while their primed version, associated to the vector field n a , are 



, = _^ + *S-^W, +Q( ^ 



and 

where 
and 



ff ' = 7 + °(^); ( 79 ) 



r = Bya (80) 



cr'o = ^(Jo ~ cro. (81) 



Let us also note that 



c'-cd-f-l (S2) 



with 



4 = ^ (83) 
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4.3 Asymptotic behavior of the Weyl components 

The leading order behavior of the Weyl components ^2, ^3 and ^f 4 satisfy the following 
relations 

^2 + ^0 ^0 + Sy^o = ^° 2 + a a + By^o, (84) 

*° =Q v a' - d vP ' 

a - _l. V x - ^ ( dvV V3vV \ - U „ ( 85 ) 

= - Oya + — O v a + I — I <7 - -OyAy 



and 



00 



2 



= -^o + 3-a + I -- 3— Ja (gg) 

V V 2 V 3 V 2 ' 

4.4 The asymptotic gauge freedom 
4.4.1 Coordinate and tetrad transformations 

Let us consider the main gauge freedom admitted in our calculation which is of the form 

fi =«(„,c,3 + *M +0 ^y (87 ) 



; = , " '. + O (r°) , 



C-C-of-l. (89) 



with d > 0. The possible further transformation of the coordinates of the sphere ((, () 
into itself, will complicate the discussion unnecessarily. 
The condition g ur = 1 imposes the relation 

w = a. (90) 

This asymptotic coordinate transformation is associated to a corresponding null tetrad 
transformation; which in the leading orders is given by 

I = du = a du + d( + d( + O ^-^ 

. 5ya _ d v a /l x 
= ct t m m + U - 



n - — + O f-^ - - — + O f- 

, <<»> 

= -n+O - 
a V r 
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^ _ V2P d Q (l\_ V2P Vw ( a c d d_ \ Q (\_ 
r d( \r 2 J r \ a du dC, 



_ - - (93) 

V2P Vwa c Vw (V 
n H — — — m + U — 



.2 



r d V \r 

since the metric expressed in terms of the new null tetrad must coincide with the metric 
expressed in terms of the original null tetrad, it is deduced that 

V V 

V=- = -; (94) 
w a 

therefore 

<5ya ^ ( 1 \ 

m = m n + U — . (95) 

ra \r 2 J 

The null tetrad transformation equations can be used to write the leading order trans- 
formation relations for the spinor dyad associated to the null tetrad [HHP73] ; namely 



o A = Va~[o A -^L A ) (96) 
ra 



" A - 1 l a . (97) 



and 

l A = - 
\/a 

Taking into account higher order transformations would include an equation of the form 

l A = ^-( L A + ho A ); (98) 
Va 

where in principle h could be of order O (r°). 

Analogous equations to (|22|) - (}24"j) relate the regular dyad at future null infinity with 
the spacetime ones, namely 

b A = n- 1 A , (99) 

i A = i A . (100) 
Then, the regular dyad at future null infinity is given by 

S A = rr 1 ~o A = ^Va-(o A -^A = ^=(d A -^ i A ) , (101) 



where we are using Q, — -. 



w \ ra I v« V a 



l A = Z A = ^- t A ; (102) 
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4.4.2 Transformation of V% and ^ 

We can now easily calculate the component ^2 of the Weyl tensor, in leading orders, with 
respect to the new null tetrad, obtaining 

*S = n-^ABOD^PtPt? = 4 (*S - 2—^3 + *2l • (103) 

cr \ a cr J 

Similarly for ty® one has 



*S = i f *? - 3^*° + 3^*° - *o ) (104) 



4.4.3 Transformation of some spin coefficients 

Recalling that the shear is defined by 



a =o A I A 'o B V AA ,o B = ^ + o(±): (100) 



then, using 

A = -^ (106) 

r d 

one obtains, for the leading orders behavior of the new shear 
a =^d (o A + A l a ) L A ' + h cr 4 ') (o B + A t B ) 

Vaa'V^(o b + Al b ) + o(J^J (107 ) 



=a (cr 



- 5 m A - A |d' A +At - A 2 // - AV) + O (J$j 



where we have used \d' to denote the thorn primed operator, k — 0, and where A is 
recognized as a quantity of {p, g}[GHP73 type {2,0}. This equation implies that 

O"0 = -T ( °0 + ^Sya + -^rd (Sytt) 2 ^rSytt 5yd (SyCt) 2 — ) (10£ 

a \ a cr cr cr k 

In order to calculate the transformation of other spin coefficients, as for example p', 
a' and r, one would need to determine the quantity h. However since the leading order 
behavior of these spin coefficients can be put in terms of V and a , one can calculate its 
transformations from those of the last two quantities. 

4.4.4 Transformation of the edth operator 

Let us recall that the edth operator can be defined in terms of the conformal factor 
P = V Pq of the sphere, by the relations 

Q v f = V2P 1 - s ^(P s f) (109) 
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and 

8 v f = V2P 1+S ^L (P~ s f) ; (110) 

where s is the spin weight of the quantity /. 

Then one can prove that the transformed edth operator acts as 



[in) 



4.4.5 Transformations to a Bondi system 

If the target coordinate system is a Bondi system, then one has 

V = l (112) 

or equivalent ly 

a = V. (113) 
With these relations the previous transformations are now 

1 ( T n 2 5Va: Tn (cVa) 2 



*S = ^(*S-^P*S + ^*JJ, (114) 

a = ^ + ^5y« - (Sya) 2 - ^Sya 3yVj , (116) 

By/ = ^vf ~ ^vV f - ^ v a f. (117) 
4.5 Rest frame and center of mass conditions 

In order to study the condition for a section to be a rest frame we need to consider the 
transformation properties of the supermomentum integrand \I/ defined in a Bondi system 
by 

= + °B0 °B0 + 9|(XB0; (118) 

where the appearance of the subindex B is used to emphasize that the quantities refer to 
a Bondi system. 

If S is a nice section then, on S one has 

^ B = —K 3 M, (119) 

where K is given by equation (jSJ . The idea is to consider the original tetrad and coordinate 
system at scri (u, C, C) as coinciding with a one parameter family of nice sections, defined 
by the conditions u = constant. 
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The expression of the supermomentum psi ^>b in terms of the intrinsic system is rather 
complicated, namely 

y B V~ 3 - 2 (Q V V) 2 a V~ 5 + Q 2 V V a V~ 4 + Q 2 v d 2 v a V~ 4 

- 4 B v d v a 8 y cVy V~ 5 - 2 d v d 2 v a d v V V' 5 - 2 Q V V B 2 v Q v a V~ 5 
+ 8 S V V B v d v a B V V V- 6 - 2 (d v V) 2 a V~ 5 - 2 d v a W V 3 V V V~ 5 

+ 8 d v a B v d v V B V V V- 6 - 8 d v a d v V {B V V) 2 V~ 7 + 2 5 y a d v V % V' 6 

+ d 2 v V do V~ 4 + d 2 v a V~ 3 - %V {d v a) 2 V~ 6 - 4 d v V B v d v a B v a V~ 6 

+ 6 5V V Q V V (B y a) 2 V' 7 + 4 d v V d v a B v a Q V V V~ 7 

+ 2 5V V Q v a (B v a) 2 V~ 7 - V a a - 2 V B v a B V V a V~ 6 

-2V (BySya) 2 V~ 6 ~2V B v d 2 v a B v a V~ 6 + 12V Q V V B v B v a S v a V~ 7 

- 10 V (d v V) 2 {p v af V- 8 + 2V Q v a Q v a V~ 4 + 8V Q v a BycVU Q v a V~ 7 

+ 4 V d v a B v d v a d v V V~ 7 -16V d v a d v V d v a B V V V~ 8 ^ 12 °^ 

- V (d v a) 2 W V V V- 7 - 2 V (Q v a) 2 Q V V Q v a V~ 8 + 2V d 2 v V (d v a) 2 V~ 7 

- 8 V B V V S v a (B y a) 2 V~ 8 + 2V 2 (B v a) 2 a V~ 7 

- 8 V 2 S v a B v S v a B v a V~ 8 + 16 V 2 Q v a S V V (S v a) 2 V~ 9 
+ V 2 (Q v a) 2 (j V- 7 + 12 V 2 (Q v a) 2 B v a Q V V V' 9 

- 10 V 3 (5ya) 2 (5ya) 2 V" 10 - V (8 y a) 2 a + 4 V 5y« SySya B v a V~ 7 

- 6 V 5Va Q V V (B v a) 2 V~ 8 - 4 V (d v a) 2 d v a S V V V~ 8 
+ 8VV (d v a) 2 (B v a) 2 V~ 9 - V (3 v a) 2 (B v a) 2 V~ 8 

+ V- 3 -cr V (cVa) 2 V~\ 

However, it is important to note that the function a which generates the transforma- 
tion (j87|) - (JH§j) is defined up to an integration constant coming from the condition ()113|) . 
Therefore at each u = uq we can choose the integration constant such that a(u = uq) = 0. 
In other words, given uq, one can always define 

a Uo (u,CC)= j U V{u\CX)du'. (121) 

Ju 

Using this choice, equation ()120|) simplifies considerably; and the nice section equation 
becomes 

2 ( 122 ) 
- — 5 ((d v V) 2 a + (d v V) 2 a ) 

= - K 3 M, 



where we are using as a the expression coming from ()121j) . 

This constitutes the condition that must be satisfied at each section of a system to 
represent a one parameter family of nice sections. 
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The center of mass condition is equation (|H7j) which in the present notation means 



Qs(w) = / C = 4 / (-WB2 + 2a B0 Q B a m + Q B (a B0 a B0 )) 

Js Js y (123) 

+ 2w m ^ B ^j dS 2 B + c.c. = 

for all a = a and b = 0; and where we are already using the nice section condition. 

All these quantities are evaluated with respect to the Bondi system, so that in terms 
of the intrinsic system one must use the transformation laws of equations (jll4|) - (jll7j) . 

5 Final comments 

In this work we have shown how to circumvent the difficulties inherent in the structure 
of the BMS group by the use of the nice section construction; in particular in the task of 
defining the intrinsic angular momentum in radiating spacetimes. We have provided, as 
far as we know, with the only definition of angular momentum which satisfies the expected 
physical properties, as those discussed at the end of section El 

In particular, let us note that detectors of gravitational waves can be considered as 
observers at future null infinity. Our construction provides for each point p at scri with a 
prescription that singles out the unique center of mass section containing p; and also with 
the appropriate rest frame to calculate intrinsic quantities like angular momentum and 
multipole moments; which are important for the description of the gravitational waves 
that one wishes to detect. 

It is worth to mention the relation of our work with previous presentations of the 
notion of angular momentum; so next we briefly comment on some (arbitrary, incomplete) 
selection of previous references. 

Early works had the ideology suggested by the conservations laws due to Komar [Kom 59j: 
where associated with a symmetry v a one can define the 2-form A ab = V[ a i>&], and define 
the 'charge' K(v) jWm80j on a 2-surface S by 

Ks(v) = ~ f A*- (124) 
8tt Js 

where A* is the dual of A. This has the same form as the expression that gives the electric 
charge q of the electromagnetic field F ab , namely 

q W = ~h \ F *- 

°n Js 

If S is a hypersurface that has as boundary S, then from Stokes' theorem one has 

K s (v) = ~ f dA*; 
8vr J s 

the fact that this holds for any such hypersurface E, is sometimes referred to as a conser- 
vation law. Noting that dA* abc = |e abcd V e A ed , one can see that if v a is a Killing symmetry 
then dA* abc = \e a bcdRtv e and therefore, assuming Einstein equations, the Komar charge 
can be related to the matter content of the spacetime. If the sources are bounded by a 



19 



world tube, then for any surface S surrounding the world tube, the Komar charge will give 
the same value. So, in the presence of a Killing symmetry, the Komar charge provides a 
simple useful notion of a conserved quantity; in particular for axis symmetric spacetimes, 
it provides with the notion of angular momentum(one component). 

Although in a general spacetime one will not have a Killing symmetry, the Wini- 
cour approach[TW66j[Win68j|Win8Qj uses the Komar charge idea, at future null infin- 
ity, associated to the asymptotic symmetries of an isolated system. When expressing 
the integral in terms of the asymptotic fields, there appears a term with the factor 
[2\&5 — 2a <Sao — <3(<7oO"o)] plus another term containing the supermomentum part. One 
must first be aware that the curvature spin coefficients definitions has been varying, and in 
some cases it is even difficult to trace the conventions; in particular the relative negative 
sign seems to be the consequence of difference with our convent ions |G HP 73 1 : however, 
one can note the non-trivial appearance of a factor of 2 in front of ty®. The Winicour 
charges were further developed in reference [GW81 where a flux was calculated. The 
relation with the symplectic fluxes was worked out in reference [AS81 . 

Using a gauge theory approach, Bramson[Bra75j, presents a notion of angular mo- 
mentum in which its integral has the factor [ty® — 2aoQa — Q(aoa )} and where no term 
associated with supermomentum appears. In fact in Bramson's approach no use of the 
generators of supertranslations is done. As we have mentioned it is essential to make use 
of this term in order to be able to pick up the notion of intrinsic angular momentum. 

Inspired by the Komar approach, Prior presented |Pri77j a definition of angular mo- 
mentum that resembles the Winicour expressions but lacking the supermomentum term. 

Using a completely different approach, Streubel |Str78j found the same expression for 
the angular momentum as Bramson's, but now with a supermomentum term. 

Penrose presented [Pen82 a notion of quasi-local mass and angular momentum asso- 
ciated with any sphere in the interior of the spacetime; making use of charge integrals of 
the Riemann tensor. His construction, when extended to a sphere at future null infinity, 
gives an integral expression where only curvature terms appear. However, as Dray and 
Streubel have shown[DS84j, one can add the shear terms with some arbitrary factors, due 
to some identities. They also presented a charge associated to generators of an unfortu- 
nately complex BMS Lie algebra. Their choice is such that the angular momentum factor 
agrees with Winicour's one. 

Recent works [KL97| |WZ00] have concentrated on the relation between angular mo- 
mentum and Lagrangian formulations, although they use background geometries in their 
constructions. 

In reference jRiz98j Rizzi presented a definition that tackles the problem of super- 
translation ambiguities. However, there is no mention in this work of the issue of super- 
momentum, and therefore it is unable to set the center of mass frame; which is needed to 
define the intrinsic angular momentum. 

Independently from the nature of the approach all these definitions have the character- 
istics that they provide with a notion of angular momentum associated to a given Lorentz 
rotation BMS generator v a and a given section S of future null infinity. The choice of a 
Lorentz rotation can be associated with a choice of a particular Bondi frame, because to 
each Bondi frame there is a Lorentz subgroup of rotation that leave the origin [u = 0) 
of the Bondi frame unchanged. Since among Bondi frames there is the supertranslation 
freedom, there is an ambiguity associated with this choice. But also any section 5" can 
be thought of as a section u = 7(C> C) °f a given Bondi frame; therefore there is also 
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a supertranslation choice associated with the election of section S. This becomes very 
important when one wants to determine the intrinsic angular momentum of the sources. 
All the definitions mentioned previously in this section suffer from this supertranslations 
ambiguities. 

In reference jMor 86j we presented a definition that fixed this problem; that required a 
preferred choice of Lorentz generators based on the properties of the spacetime in the limit 
for u — > — oo; that is in the limit to spacelike infinity. Therefore given any S we provided 
with a charge integral of the Riemann tensor that had the correct angular momentum 
behavior, with an immediate flux law and without supertranslation ambiguities. The 
remaining difficulty with our previous work was that it would not necessarily pick up the 
natural expected notion of intrinsic angular momentum of the sources for the cases of 
radiating spacetimes. 

From the naive comparison of the mathematical expression for angular momentum 
factors on a given section 5*, one could distinguish mainly the following cases: 



Angular momentum factor 


Supermomentum factor 


appearing in references 


Winicour type 
Winicour type 
Winicour type 
Winicour type 
Bramson type 
Bramson type 


Winicour supermomentum 
absent 

Geroch supermomentum 
psi supermomentum 
absent 

psi supermomentum 


[TW66j[Win68 Win80j 
|Pr777] 

p^jwznn] 

|Bra?5| |H,iz98j 

jStr78j lMor86j [present] 



Actually Rizzi's integrand does not coincide with Bramson's exactly, but it is closer to this 
than to Winicour type integrands. The fact that two or more references appear in the same 
row, does not mean that they are equivalent. This table only shows the similarities of the 
factors involving the curvature tensor; but the definitions are completed with the factors 
involving the information of the BMS generator that one is using. In particular, references 
|Str78j and [Mor86j both are inspired from Penrose work |Pen82j : but applying the defini- 
tion of reference |Str78j to two sections Si and S2, due to the supertranslation ambiguity 
problem, it is not clear how to relate the two quantities. Instead reference [Mor86j consti- 
tutes an extension of Penrose expression to the whole of scri, making sense for radiating 
spacetimes, solving the ambiguities of supertranslations by a physical condition and pro- 
viding with an immediate flux law, which is zero in the absence of gravitational radiation. 
It should be stressed that the Penrose expressions can be made to agree to Winicour or 
Bramson type angular momentum factors, equivalently to agree with Geroch or psi super- 
momentum factors. It seems that in reference |Str78j the choice was intended to get the 
Bramson angular momentum type, although later the author changed his mind for the 
Winicour type |DS84j : instead in reference [Mor86j the stress was on the supermomentum 
part, due to the properties of the psi supermomentum[Mor86j[Mor88j[MD98 j DM00 , and 
as a consequence the angular momentum factor was obtained to be of Bramson type. 

Our present construction; which can also be considered a development of the Penrose 
approach, and a mayor improvement over reference |Mor86j . provides for the first time 



with a global well behaved notion of intrinsic angular momentum for radiating sources, 
and also gives the notion of instantaneous center of mass frame; which is essential for the 
discussion of higher multipole moments. 

For example, if one is studying a system whose gravitational radiation is appropriately 
described by the quadrupole radiation formula, at an instant of time; then in a long lapse 
of time, the corresponding center of mass frame will be supertranslated with respect to 
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the original one. Therefore, in order to know what is the quadrupole moment at a later 
time, one needs to settle what is the center of mass frame first; which we have done here. 

We next review some expected properties that one may require a charge integral of 
BMS symmetries v a to satisfy. The first [3 are adapted from reference [AW82J and we 
have added another one. 

1. Qs(v) should be linear in the BMS symmetry v a . 

2. The expression Qs{v) should only involve local fields at 5* on scri. 

3. For the case when v a is a translation, Qs{v) should agree with the corresponding 
component of the Bondi momentum evaluated at S. 

4. If v a is the natural extension to scri of a Killing symmetry of the spacetime then 
Qs(v) should be independent of S, if the vacuum Einstein equations are satisfied en 
a neighborhood of scri, and it should be proportional to the Komar charge Ks{v) 
of equation (JT2H). 

5. In Minkowski space, Qs{v) should vanish for all BMS symmetries v a and section S. 

6. Given a symmetry v a there should exist a local flux F(v), linear in v a , such that 
Qs{v) — Qs'(v) = J s F(v); where S is the region on scri bounded by S and S'. 

7. The flux F(v) should vanish in absence of gravitational radiation. 

8. Having three consecutive stages at scri, a first stationary region Si, a second radiat- 
ing stage S2 and a third stationary region S3, then there is a continuous prescription 
for the center of mass sections S cm provided by the angular momentum calculated 
from Qs{v) such that J£P b = 0, where J ab is the angular momentum calculated 
from Qs cm ( v cm), in which v cm are the Lorentz rotations leaving S cm invariant, and P b 
is the Bondi momentum. 

One can readily see that our definition of charge integrals of the Riemann tensor at scri, 
which uses the charge Qs{v) of equation (pD|). satisfies all condition except HJ although it 
is still an open question whether condition 0] is satisfied or not. Ours is the only definition 
we know to satisfy property |H1 

We have seen that the rest frame and center of mass condition for the intrinsic systems 
are rather complicated. However one should have in mind that for astrophysical interesting 
systems, the amount of gravitational radiation is small, and therefore several quantities 
can be approximated around stationary values. 

This is precisely the key idea we plan to implement in future works. 
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